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G^ i One of the most important problems in astrophysics concerns the nature of the dark 

, matter in galactic halos, whose presence is implied mainly by the observed flat rotation 

curves in spiral galaxies. In the framework of a baryonic scenario the most plausible can- 
didates are brown dwarfs, M-dwarfs or white dwarfs and cold molecular clouds (mainly 
of H2)- The former can be detected with the ongoing microlensing experiments, which 
are rapidly leading to important new results. The French collaboration EROS and the 
American-Australian collaboration MACHO have reported until August 1997 the obser- 
^v^j i vation of ~ 16 microlensing events by monitoring during several years the brightness 

^v^j ' of millions of stars in the Large Magellanic Cloud and one event towards the Small 
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Magellanic Cloud. In particular, the MACHO team found 8 microlensing candidates by 
analysing their first 2 years of observations. This implies that the halo dark matter frac- 
tion in form of MACHOs (Massive Astrophysical Compact Halo Objects) is of the order 
of 45-50% assuming a standard spherical halo model. More than 150 microlensing events 
have been detected in the direction of the galactic bulge by the MACHO, OGLE and 
DUO teams. The measured optical depth implies the presence of a bar in the galactic 
centre. Here, we give an overview of microlensing and the main results achieved so far. 
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One of the most important problems in astrophysics concerns the nature of the dark 
matter in galactic halos- .whose presence is implied by the observed flat rotation 



1 Introduction 



curves in spiral galaxies BB, the X-ray diffuse emission in elliptical galaxies as well 
, as by the dynamics of galaxy clusters. Primordial nucleosynthesis entails that most 

of the baryonic matter in the Universe is nonluminous, and such an amount of dark 
matter falls suspiciously close to that required by the rotation curves. Surely, the 
standard model of elementary particle forces can hardly be viewed as the ultimate 
theory and all the attempts in that direction invariably call for new particles. Hence, 
jJJ ■ the idea of nonbaryonic dark matter naturally enters the realm of cosmolgy and 

may help in the understanding of the process of galaxy formation and clustering of 
galaxies. _ 
The problem of .dark matter started already with the pioneering work of OortB 
in 1932 and ZwickyQ in 1933 and its mistery is still not solved. Actually, there are 
several dark matter problems on different scales ranging from the solar neighbour- 
hood, galactic halos, cluster of galaxies to cosmological scales. Dark matter is also 
needed to understand the formation of large scale structures in the universe. Many 
candidates have been proposed, either baryonic or not, to explain dark matter. 

Here, we discuss the dark matter problem in the halo of our Galaxy in con- 
nection with microlensing searches. We present the basics of microlensing and an 
overview of the results obtained so far, without, however, being exhaustive and this 
also for the references we quote. The content is as follows: first, we review the 
evidence for dark matter in the halo of our Galaxy. In Section 3 we present the 
baryonic candidates for dark matter and in Section 4 we give an overview of the 
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results of microlensing searches achieved so far. In Section 5 we discuss the basics 
of microlensing (optical depth, microlensing rates, etc.) and in Section 6 we briefly 
present a scenario in which part of the dark matter is in the form of cold molecular 
clouds (mainly of H 2 ). 

2 Mass of the Milky Way 

The best evidence for dark matter in galaxies comes from the rotation curves of 
spirals. Measurements of the rotation velocity v ro t of stars up to the visible edge 
of the spiral galaxies and of HI gas in the disk beyond the optical radius (by 
measuring the Doppler shift in the 21-cm line) imply that v ro t ~ constant out to 
very large distances, -rather than to show a Keplerian falloff. These observations 
started around 1970B, thanks to the improved sensitivity in both optical and 21- 
cm bands. By now there are observations for over thousand spiral galaxies with 
reliable rotation curves out to large radiiH. In almost all of them the rotation curve 
is flat or slowly rising out to the last measured point. Very few galaxies show falling 
rotation curves and those that do either fall less rapidly than Keplerian have nearby 
companions that may perturb the velocity field or have large spheroids that may 
increase the rotation velocity near the centre. 

There are also measurements of the rotation velocity for our Galaxy. However, 
these observations turn out to be rather difficult, and the rotation curve has been 
measured only up to a distance of about 20 kpc. Without any doubt our own galaxy 
has a typical flat rotation curve. A fact this, which implies that it is possible to 
search directly for dark matter characteristic of spiral galaxies in our own Milky 
Way. 

In oder to infer the total mass one can also study the proper motion of the 
Magellanic Clouds and of other satellites of our Galaxy. Recent studies uBu do not 
yet allow an accurate determination of v rot (LMC)/vo (v = 210±10 km/s being the 
local rotational velocity). Lin et al. 13 analyzed the proper motion observations and 
concluded that within 100 kpc the Galactic halo has a mass ~ 5.5 ± 1 x lO 11 !!© and 
a substantial fraction ~ 50% of this mass is distributed beyond the present distance 
of the Magellanic Clouds of about 50 kpc. Beyond 100 kpc the mass may continue 
to increase to ~ 10 12 Mq within its tidal radius of about 300 kpc. This value £pr the 
total mass of the Galaxy is in agreement with the results of Zaritsky et al. 0, who 
found a total mass in the range 9.3 to 12.5 xl0 n Mo, the former value by assuming 
radial satellite orbits whereas Jhe latter by assuming isotropic satellite orbits. 

The results of Lin et al. U suggest that the mass of the halo dark matter up 
to the Large Magellanic Cloud (LMC) is roughly half of the value one gets for 
the standard halo model (with flat rotation curve up to the LMC and spherical 
shape), implying thus the same reduction for the number of expected microlensing 
events. KochanekEI analysed the global mass distribution of the Galaxy adopting a 
Jaffe model, whose parameters are determined using the observations on the proper 
motion of the satellites of the Galaxy, the Local Group timing, constraint and the 
ellipticity of the M31 orbit. With these observations Kochaneklj concludes that the 
mass inside 50 kpc is 5.4 ± 1.3 x 10 11 M Q . This value becomes, however, slightly 
smaller when using only the satellite observations and the disk rotation constraint, 
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in this case the median mass interior to 50 kpc is in the interval 3.3 to 6.1 (4.2 to 
6.8) without (with) Leo I satellite in units of 10 u Mq. The lower bound without 
Leo I is 65% of the mass expected assuming a flat rotation curve up to the LMC. 

3 Baryonic dark matter candidates 

Before discussing the baryonic dark matter we would like to mention that another 
class of candidates which is seriously taken into consideration is the so-called cold 
dark matterp-which consists for instance of axions or supersymmetric particles like 
neutralinos L9. Here, we will not discuss cold dark matter in detail. However, 
recent studies seem to point out that there is a discrepancy between the calculated 
(through N-body simulations) rotation curve, for .dwarf galaxies assuming an halo 
of cold dark matter and the measured curves liilOH. If this fact is confirmed, this 
would exclude cold dark matter as a major constituent of the halo of dwarf galaxies 
and possibly also of spiral galaxies. 

From the Big Bang nucleosynthesis model EJll3 and from the observed abun- 
dances of primordial elements one infers: 0.010 < HqCIb < 0.016 or with ho ~ 0.4—1 
one gets 0.01 < Hb < 0.10 (where f2_B = PbI ' Per it, and p cr u = 3H 2 /8itG). Since for 
the amount of luminous baryons one finds f^« m <C £Ib , it follows that an important 
fraction of the baryons are dark. In fact the dark baryons may well make up the 
entire dark halo matter. 

The halo dark matter cannot be in the form of hot ionized hydrogen gas other- 
wise there would be a large X-ray flux, for which there are stringent upper limitsE3. 
The abundance of neutral hydrogen gas is inferred from the 21-cm measurements, 
which show that its contribution is small. Another possibility is that the hydro- 
gen gas is in molecular form clumped into cold clouds, as we will briefly discuss in 
Section 6. Baryons could_otherwise have been processed in stellar remnants (for 
a detailed discussion seeEO). If their mass is below ~ 0.08 M Q they are too light 
to ignite hydrogen burning reactions. The possible origin of such brown dwarfs or 
Jupiter like bodies (called also MACHOs) , by fragmentation or by some other mech- 
anism, is at present not well understood. It has also been pointed out that the mass 
distribution of the MACHOs, normalized to the dark halo mass density, cpuld be 
a smooth continuation of the known initial mass function of ordinary starslla. The 
ambient radiation, or their own body heat, would make sufficiently small objects of 
H and He evaporate rapidly. The condition that the rate of evaporation of such a 
hydrogenoid sphere be insufficient to halve its mass in a billion years leads to the 
following lower limit on their massES M > 10" 7 A/ Q (T s /30 K) 3 / 2 (l g cm^/p) 1 / 2 
(Tg being their surface temperature and p their average density, which we expect 
to be of the order ~ 1 g cto~ 3 ). 

Otherwise, MACHOs might be either M-dwarfs or else white dwarfs. As a 
matter of fact, a deeper analysis shows that the M-dwarf option looks problematic. 
The null result of several searches for low-mass stars both in the disk and in the 
halo of our Galaxy suggests that the halo cannot be mostly in the form of hydrogen 
burning main sequence M-dwarfs. Optical imaging of high-latitude fields taken with 
the Wide Field Camera of the Hubble Space Telescope indicates that less than ~ 6% 
of the halo can be in this formlia. Observe, however, that these results are derived 
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under the assumption of a smooth spatial distribution of Mj-dwarfs, and become 
considerably less severe in the case of a clumpy distributionEJ'Eil. 

A scenario with white dwarfs as a major constituent of the galactic halo dark 
matter has been explored E3. However, it requires a rather ad hoc initial mass 
function sharply peaked around 2-6 M©. Future Hubble deep field exposures 
£ould either find the white dwarfs or put constraints on their fraction in the halo 
t3. Also a substantial component of neutron stars and black holes with mass higher 
than ~ 1 Mq is excluded, for otherwise they would lead to an overproduction of 
heavy elements relative to the observed abundances. 



4 Present status of microlensing research 

It has been pointed out by PaczynskiEl that microlensing allows the detection of 
MACHOs located in the galactic halo in the mass raagetJ 10~ 7 < M/Mq < 1, as 
well as MACHOs in the disk or bulge of our Galaxyl2§Ea. Since this first proposal 
microlensing searches have turned very quickly into reality and in about a decade 
they have become an important tool for astrophysical investigations. Microlensing 
is very promising for the search of planets around other stars in our Galaxy and 
generates also very large databases for variable stars, a field which has already 
benefitted a lot. Because of to the increase of observations, due also to the fact 
that new experiments are becoming operative, the situation is changing rapidly 
and, therefore, the present results should be considered as preliminary. Within few 
years the amount of data will be such that several problems will be solved or at 
least allow to achieve substantial progress. The following presentation is also not 
exhaustive with respect to all what has been found so far. 



4.1 Towards the LMC and the SMC 

In September 1993 the French collaboration EROS0 announced the discovery of 
2 microlensing_candidates and the American-Australian collaboration MACHO of 
one candidateEa by monitoring stars in the LMC. 

In the meantime the MACHO team reported the observation of altogether 8 
events (one is a binary lensing event) analysing ±heir first two years of data by 
monitoring about 8.5 million of stars in the LMCE£l. The inferred optical depth is 
r = 2.1±o"7 x 10~ 7 when considering 6 events |] (see Table 3) (or r = 2.9±q;| x 1CT 7 
when considering all the 8 detected events). Correspondingly, this implies that 
about 45% (50% respectively) of the halo dark matter is in form of MACHOs and 
they find an average mass 0.5^q'2-M© assuming a standard spherical halo model. It 
may well be that there is also a contribution of events due to MACHOs located in 
the LMC itself or in a_thich disk of our galaxy, the corresponding optical depth is 
estimated to be only a r = 5.4 x 10~ 8 . Other events have been detected towards 
the LMC by the MACHO group, which have been put on their list of alert events. 
In particular 2 events are reported in 1996 and already 4 events in 1997. The full 
analysis of the 1996 and 1997 seasons is still not published. 

"In fact, the two disregarded events are a binary lensing and one which is rated as marginal. 
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EROS has also searched for very-low mass MACHOs hy_looking for microlensing 
events with time scales ranging from 30 minutes to 7 days El The lack of candidates 
in this range places significant constraints on any model for the halo that relies on 
objects in the range 5 x 1CP 8 < M/Mq < 2 x 10~ 2 . Indeed, such objects may 
make up at most 20% of the halo dark matter (in the range between 5 x 10~ 7 < 
M /M < 2 x 1CIZ? at most 10%). Similar conclusions have_also been reached by the 
MACHO group Ej. Recently, the MACHO team reported £3 the first discovery of a 
microlensing event towards the Small Mgellanic Cloud (SMC). The full analysis of 
the four years data on the SMC is still underway, so that more candidates may be 
found in the near future. A rough estimate of the optical depth leads to about the 
same value as found towards the LMC. 

The EROS group has completed his first campaign based on photographic plates 
and CCD, the latter one for the short duration events. Since the middle of 1996 
the EROS group has put into operation a new 1 meter telescope, located in La Silla 
(Chile), and which is fully dedicated to microlensing searches using CCD cameras. 
The improved experiment is called EROS II. 



4-2 Towards the galactic centre 

Towards the galactic bulge the Polish- American team OGLeU announced his first 
event also in September 1993. Since then OGLE found in their data from the 1992 - 
1995 observing seasons altogether 18 microlensing events (one being a binary lens). 
Based on theix first 9 events the OGLE team estimated the optical depth towards 
the bulge asLJ r = (3.3 ± 1.2) x lOz^pIThis has to be compared with the theoretical 
calculations which lead to a valueOo r ~ (1 — 1.5) x 10 -6 , which does, however, 
not take into account the contribution of lenses in the bulge itself, which might 
well explain the discrepancy. In fact, when taking into account-also the effect of 
microlensing by galactic bulge stars the optical depth gets biggerEa and might easily 
be compatible with the measured value. This implies the presence of a bar in the 
galactic centre. In the meantime the OGLE group got a new dedicated 1.3 meter 
telescope located at the Las Campanas Observatory. The OGLE-2 collaboration 
has started the observations in 1996 and is monitoring the bulge, the LMC and the 
SMC as well. 

The French DUO Ell team found 12 microlensing events (one of which being a 
binary event) by monitoring the galactic bulge during the 1994 season with the ESO 
1 meter Schmidt telescope. The photographic plates were taken in two different 
colors to test achromaticity. The MACHO Ell collaboration found by now more 
than ~ 150 microlensing events towards the galactic bulge, most of which are listed 
among the alert events, which are constantly updated |. They found also 3 events 
by monitoring the spiral arms in the region of Gamma Scutum. During their first 
season they found 45 events towards the bulge. The MACHO team detected also 
in a long-duration event the parallax effect due to the motion of the Earth around 
the SunEZL The MACHO first year data leads to an estimated optical depth of 
r ~ 2.43^q'45 x 10~ 6 j which is roughly in agreement with the OGLE result, and 



^ Current information on the MACHO Collaboration's Alert events is maintained at the WWW 



site http:/ /darkstar. astro. washington.edu 
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which also implies the presence of a bar in the galactic centre. These results are 
very important in order to study the structure of our Galaxy. In this respect the 
measurement towards the spiral arms will give important new information. 

4-3 Towards the Andromeda galaxy 

Microlensing searches have also been conducted towards M31. In this case, however, 
one has to use the so-called "pixel-lensing" method, since the source stars are in gen- 
ial no longer resolvable. Two groups have performed searches: the French AGAPE 
E using the 2 meter telescope at Pic du Midi and the American VATT/COLUMBIA 
23, which used the 1.8 meter VATT-tclcscope located on Mt. Graham and the 4 
meter KNPO telescope. Both teams showed that the pixel-lensing method works, 
however, the small amount of observations done so far does not allow to draw firm 
conclusions. The VATT/COLUMBIA team found six candidates which are con- 
sistent with microlending, however, additional observations are needed to confirm 
this. Pixel-lensing could also lead to the discovery of microlensing events towards 
e M87 galaxy, in which case the best would be to use the Hubble Space Telescope 
. It might also be interesting to look towards dwarf galaxies of the local group. 

4-4 Further developments 

A new collaboration between New Zealand and Japan, called MOA, started in june 
1996 to-perform observations using the 0.6 meter telescope of the Mt. John Obser- 
vatory ell. The targets are the LMC and the galactic bulge. They will in particular 
search for short timescale (~ 1 hour) events, and would then be particularly sensitive 
to objects with a mass tipical for brown dwarfs. 

It has to mentioned that there are also collaborations between different obser- 
vatories (for instance PLANETc3 and GMANE3) with the aim to perform accurate 
photometry on alert microlensing events. The GMAN collaboration was able to 
accurately get photometric data on a 1995 event towards the galactic bulge—The 
light curve shows clearly a deviation due to the extension of the source starcj. A 
major goal of the RLAMET and GMAN collaborations is to find planets in binary 
microlensing eventsc§c3o. Moreover, microlensing searches are also very powerful 
ways to get large database for the study and discovery of many variable stars. 

At present the only information available from a microlensing event is the time 
scale, which depends on three parameters: distance, transverse velocity and mass 
of the MACHO. A possible way to get more information is to observe an event from 
different locations, with typically an Astronomical Unit in_ser>aration. This could 
be achieved by putting a parallax satellite into solar orbit c3o. 

The above list of presently active collaborations and main results shows clearly 
that this field is just at the beginning and that many interesting results will come 
in the near future. 

5 Basics of microlensing 

In the following we present the main features jq£ -microlensing, in particular its 
probability and rate of events (for reviews see alsoE30, whereas for double lenses see 
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for instance ref. An important issue is the determination from the observations 
of the mass of the MACHOs that acted as gravitational lenses as well as the fraction 
of halo dark matter they make up. The most appropriate way to compute the 
average mass and other important information is to use the method of mass moments 
developed by De Rujula et al. 13, which will be briefly discussed in Section 5.6. 



5.1 Microlensing probability 

When a MACHO of mass M is sufficiently close to the line of sight between us 
and a more distant star, the light from the source suffers a gravitational deflection. 
The deflection angle is usually so small that we do not see two images but rather a 
magnification of the original star brightness. This magnification, at its maximum, 
is given by 

_ u* + 2 

max ~ U (w 2 +4)V2 ' [L > 

Here u = d/ Re (d is the distance of the MACHO from the line of sight) and the 
Einstein radius Re is defined as 

AGMD . . 

R%= —x(l-x) (2) 

c z 

with x = s/D, and where D and s are the distance between the source, respectively 
the MACHO and the observer. 

An important quantity is the optical depth r opt to gravitational microlensing 
defined as 

T op t = I dx — —p(x)D 2 x(l - x) (3) 
Jo c 

with p(x) the mass density of microlensing matter at distance s = xD from us 
along the line of sight. The quantity r op t is the probability that a source is found 
within a radius Re of some MACHO and thus has a magnification that is larger 
than A = 1.34 (d < R E ). 

We calculate T opt for a galactic mass distribution of the form 

P(r) = Ma2 2 t? c) . (4) 

| v | being the distance from the Earth. Here, a is the core radius, po the local dark 
mass density in the solar system and Roc the distance between the observer and 
the Galactic centre. Standard values for the parameters are po = 0.3 GeV/cm 3 = 
7.9 10~ 3 Mq/pc 3 , a — 5.6 kpc and Rgc = 8.5 kpc. With these values we get, far a 
spherical halo, T opt ~ 5 x 1CT 7 for the LMC and r opt ~ 7 x lfr 7 for the SMCEl 

The magnification of the brightness of a star by a MACHO is a time-dependent 
effect. For a source that can be considered as pointlike (this is the case if the 
projected star radius at the MACHO distance is much less than Re) the light curve 
as a function of time is obtained by inserting 

He 
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into eq.(y), where vt is the transverse velocity of the MACHO, which can be inferred 
from the measured rotation curve (vt ~ 200 km/ s). The achromaticity, symmetry 
and uniqueness of the signal are distinctive features that allow to discriminate a 
microlensing event from background events such as variable stars. 

The behaviour of the magnification with time, A(t), determines two observables 
namely, the magnification at the peak A(0) - denoted by A max - and the width of 
the signal T (defined as being T — Re/vt)- 



5.2 Microlensing rate towards the LMC 

The microlensing rate depends on the mass and velocity distribution of MACHOs. 
The mass density at a distance s — xD from the observer is given by eq.(^). The 
isothermal spherical halo model does not determine the MACHO number density 
as a function of mass. A simplifying assumption is to let the mass distribution 
be independent of the position in the galactic halo, i.e., we assume the following 
factorized form for the number density per unit mass dn/dM, 

dn dno a 2 + Rr<r dug 

dM dM = -dHa 2 + Rl c + D^-2DR GcX co S a dfi = ' (6) 

with p = M/Mq (a is the angle of the line of sight with the direction of the 
galactic centre, which is 82° for the LMC), uq not depending on x and is subject to 
the normalization J dp^-M — po- Nothing a priori is known on the distribution 
dno/dM. 

A different situation arises for the velocity distribution in the isothermal spher- 
ical halo model, its projection in the plane perpendicular to the line of sight leads 
to the following distribution in the transverse velocity vt 

f(v T ) = 4-VTe- v * /v ». (7) 



v 



(vh ~ 210 km/ s is the observed velocity dispersion in the halo). 

In order to find the rate at which a single star is microlensed with magnification 
A > A m i n , we consider MACHOs with masses between M and M + SM, located 
at a distance from the observer between s and s + Ss and with transverse velocity 
between vt and vt + Svt- The collision time can be calculated using the well-known 
fact that the inverse of the collision time is the product of the MACHO number 
density, the microlensing cross-section and the velocity. The rate dT, taken also as 
a differential with respect to the variable u, at which a single star is microlensed in 
the interval dLidudvTdx is given byrnri 

dT = 2v T f(vT)Dr E \Lix(l - x)} 1/2 H(x)^-dududv T dx, (8) 

dp 

with 

r% = „ ( 3. 2 x 10 9 fcm) 2. (9) 

& 

One has to integrate the differential number of microlensing events, dN ev = 
N^tobsdT, over an appropriate range for fx, x, u and Vt, in order to obtain the total 
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number of microlensing events which can be compared with an experiment moni- 
toring N+ stars during an observation time i t, s and which is able to detect a mag- 
nification such that A max > Ath- The limits of the u integration are determined 
by the experimental threshold in magnitude shift, Atuth' we have < u < uth- 

The range of integration for /i is where the mass distribution dno/d/i is not 
vanishing and that for x is < x < Dh/D where Dh is the extent of the galactic 
halo along the line of sight (in the case of the LMC, the star is inside the galactic halo 
and thus Dh/D = 1.) The galactic velocity distribution is cut at the escape velocity 
v e sa 640 km/ s and therefore vt ranges over < vt < v e . In order to simplify the 
integration we integrate vt over all the positive axis, due to the exponential factor 
in f(vx) the so committed error is negligible. 

However, the integration range of dfidudvxdx does not span all the interval we 
have just described. Indeed, each experiment has time thresholds T m i n and T max 
and only detects events with: T m i n < T < T rnax , and thus the integration range 
has to be such that T lies in this interval. The total number of micro-lensing events 
is then given by 

N ev = J dN ev e(T) , (10) 

where the integration is over the full range of dfidudvTdx. e(T) is determined ex- 
perimentally Ej'EJ or if not known can, for instance, be taken as follows e(T) ~ 
Q(T — T m i n )Q(T ma x — T). T is related in a complicated way to the integration vari- 
ables, because of this, no direct analytical integration in eq.([Io|) can be performed. 
To evaluate eq. (|l0|) we define an efficiency function e (fi) 

,_ JdN* v (fl) e(T) 
J dN* v (n) 

which measures the fraction of the total number of microlensing events that meet 
the condition on T at a fixed MACHO mass M — /2Mq. We now can write the 
total number of events in eq. ( |l0| ) as 



dN ev e (/i) • (12) 

Due to the fact that e is a function of /i alone, the integration in d^dudvTdx 
factorizes into four integrals with independent integration limits. 
The average lensing duration can be defined as follows 

<T>=i J dTT(x,fi,v T ) , (13) 

where T(x, n, vt) = Re{x, fi)/vT- One easily finds that < T > satisfies the following 
relation 

<T>=% u TH . (14) 

7Tl 

In order to quantify the expected number of events it is convenient to take as 
an example a delta function distribution for the mass. The rate of microlensing 
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events with A > A m i n (or u < u max ), is then 

1 f 1 

r(j4 mi „) = u maa f = u max Dr E vhyt-—^ / rfx[a;(l - a;)] 1/2 iZ"(x) . (15) 

M VA 1 Jo 



Inserting the numerical values for the LMC (D=50 kpc and a = 82°) we get 

vh \ I Po \ 1 



T = 4 x 10 



-13 



210 km/s ) V0.3 GeV/cm 3 J y/M/M Q 



(16) 



For an experiment monitoring TV* stars during an observation time t b s the total 
number of events with a magnification A > A m i n is: N ev (A m i n ) — N- k t f )s T{A m i n ). 
In the following Table 1 we show some values of N ev for the LMC, taking t h s = 1 
year, TV* = 10 6 stars and A min = 1.34 (or Am min = 0.32). 

Table 1 



MACHO mass in M 


Mean Re in km 


Mean microlensing time 




lO" 1 


0.3 x 10 9 


1 month 


4.5 


lO" 2 


10 8 


9 days 


15 


io~ 4 


10 7 


1 day 


165 


10~ 6 


10 6 


2 hours 


1662 



As mentioned in Sect. 4.1 the MACHO team found (till September 1997) 
altogether 14 events (one of which being a binary lens) towards the LMC and one 
event towards the SMC. The EROS team found 2 events using photographic plate 
techniques. 

5.3 Microlensing rate towards M31 

Gravitational mkxolensing is also useful for detecting MACHOs in the halo of 
nearby galaxiesc3i53 such as M31 or M33, for which the experiments AGAPE and 
VATT/Columbia are under way. In fact, it turns out that the massive dark halo 
of M31 has an optical depth to microlpasmg .which is of about the same order of 
magnitude as that of our own galaxy Ej'EJjEa ~ 10~ 6 . Moreover, an experiment 
monitoring stars in M31 is sensitive to both MACHQa in our halo and in the one 
of M31. One can also compute the microlensing rateEa for MACHOs in the halo of 
M31, for which we get 



T = 1.8 x 10" 



vh 



P(0) 



1 



(17) 



210 km/s J \1 Gev/cm 3 J ^/M/Mq 

(p(0) is the central density of dark matter.) The average lensing time is given by 

< T >~ (125 days) y/M/M e . (18) 

In the following Table 2 we show some values of N£ v due to MACHOs in the 
halo of M31 with t b s = 1 year and TV* = 10 6 stars. In the last column we give the 
corresponding number of events, N evi due to MACHOs in our own halo. The mean 
microlensing time is about the same for both types of events. 
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Table 2 



MACHO mass in M 


Mean Re in km 


Mean microlensing time 


N a 

ev 


N ev 


lO" 1 


7 x 10 s 


38 days 


6 


3 


1(T 2 


2 x 10 s 


12 days 


21 


12 


lO" 4 


2 x 10 7 


30 hours 


210 


129 


10~ 6 


2 x 10 6 


3 hours 


2100 


1290 



N£ v is almost by a factor of two bigger than N ev (see alsoE3). Of course these 
numbers should be taken as an estimate, since they depend on the details of the 
model one adopts for the distribution of the dark matter in the halo. To distinguish 
between events due to a MACHO in our halo or in the one of M31 might be rather 
difficult. 



5.4 Microlensing rate towards the galactic bulge 

We compute the microlensing rate T for an experiment monitoring stars in the 
galactic bulge in Baade's window of galactic coordinates (longitude and latitude): 
I = 1°, b = — 3.°9. D = 8.5 kpc is the distance to stars in the galactic bulge, 
te = 1.25 x 10 9 km and we use the definition T = Re/vt- ■—■ 
The number density of disk stars per unit mass is given byEj 

dn dno ( Dx I sin b I Dx cos b\ driQ . . 

dM = dM GXP {- 300 pc + T5kp^J = dM Hd[x) ' (19) 

where the galactic longitude I = 0° has been adopted and = w ^fi 
Pd(M) = O.O5M pc -3 ; Hd can be written as follows 

fxD (1 - 11.7 I sin b |)\ 
g ^ =eX H 3.5 kpc -)> (20) 

using also the fact that cos b w 1 for the galactic bulge. 

For the contribution H (x) of the halo dark matter in form of MACHOs located 
in the disk one uses the distribution for the number density per unit mass dn/ dM as 
given in eq. (||) with a ~ 4° (the angle between the line of sight and the direction 
of the galactic centre). 

With the above distribution for the lenses in the disk it follows that the corre- 
sponding optical depth is r ~ 7x 10 -7 (and r ~ 1.2x 10~ 7 for the halo contribution). 

In computing V one must also take into account the fact that both the source 
and the observer are in motionE3. Relevant are only the velocities transverse to the 
line of sight. The transverse velocity of the microlensing tube at position xD is: 
vt(x) = (1 — x)vqj_ + xv s ±, and its magnitude is 

v t (x) = y/(l - xf I v Q± I 2 +x 2 I v s± I 2 +2.t(1 -x) I v Q± \\ v s± \ cos 6 , (21) 

where v s _i and vq±_ are the source and the solar velocities transverse to the line of 
sight and 8 the angle between them. 
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For the velocity distribution of the MACHOs or the faint disk stars we consider 
an isothermal spherical model, which in the rest frame of the galaxy is given by 

f(v)d 3 v = ^— e-* 2 l* 2 H d 3 v . (22) 

Since only the transverse velocities are of relevance cylindrical coordinates can be 
used and the integration made over the velocity component parallel to the line of 
sight. Moreover, due to the velocities of the observer and the source, the value of 
the transverse velocity gets shifted by Vt(x). The distribution for the transverse 
velocity is thus 

f(v T )dv T = -Lr- e-Wr-*) 8 /** VT dv T , (23) 

where vh ~ 30 km s _1 is the velocity dispersionEl 

The random velocity of the source stars in the bulge are again described by an 
isothermal spherical distribution, whose transverse velocity distribution is 

g(v sl _)dv s± = — =- e - v2 ^ /v v v s± dv s± , (24) 

where vp = 156 km s^ 1 is the velocity dispersion 

Taking all the above facts into account, T turns out to bel£j 

(f°° dn \ f 27T f°° 
J ^ ~djjL d ^J J dd J dVsl ~ 9 ^ Vs± ^ 



y/x(l-x) Hi(x) e-* / dy y 2 I (2yr)) e~ v , (25) 



where y = rj(x, 9, v s ± , Vq±) = and Iq is the modified Bessel function of 
order 0. In the limit of stationary observer and source star (vq± = v s ± = 0), 
r\ = and Iq = 1; Hi means either Hd or H; uth is related to the minimal 
experimentally detectable magnification Ath = A[u = uth}- vq± is ~ 220 km s _1 
(more precisely it should be multiplied by cos I, where I is the galactic longitude 
but since 1 = 1°, cos I w 1). In computing T one should also take into account 
the limited measurable range for the event duration T, which translates into a 
modification of the integration limits. A fact that can be described by introducing 
an efficiency functioned eo(n). 

Assuming a delta-function-type distribution for the masses 

dn _ p 6(n-fl) 
dfi Mq h 

eq. ( p5| ) can be integrated. With = 10 6 stars and t bs = 1 year one gets 

yf\x \5 x 1O _2 M pc^- 3 / 



for Hi = Ha, and 



N ev = ^r ( if) UTH, (2M 
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for Hi = H. The numerical factor in eq. ( p7j ) for N ev as a function of b, the galactic 
latitude, varies between 6.8 for 6 = 0° and 1.9 for b = 5°, whereas the factor for 
N ev of eq. fl2q) remains practically unchanged. 

As mentioned in Section 4.2 the results clearly show that there is a bar in 
the galactic centre and that one has to consider also the bar-bar (or bulge-bulge) 
contribution in order to explain the observations. 



5.5 Most probable mass for a single event 

The probability P that a microlensing event of duration T and maximum ampli- 
cation A max be produced by a MACHO of mass /i (in units of M Q ) is given by 

P(M, T) cx £ jf 1 rix(x(l - x)) 2 ff (x)^ (- X) ) , (29) 

which does not dependent on A maa; and P(/U, T) = P(/j,/T 2 ). The measured values 
for T towards the LMC are listed in Table 3, where [imp is the most probable value. 
The normalization is arbitrarily chosen such that the maximum of Pt ' n 'mp , T) = 
1. We find that the maximum corresponds to fir% / 'v 2 H T 2 = 13.0 @fl The 50% 
confidence interval embraces for the mass fi approximately the range IIZ^mp up 
to 3fiMP- Similarly one can compute P(/x, T) also for the bulge eventsEJ. 



Table 3: Values of ^mp (in Mq) for eight microlensing events detected in the 
LMC (Ai = American-Australian collaboration events (i = 1,..,6); F\ and F2 French 
collaboration events). For the LMC: vh — 210 km s _1 and rE = 3.17 x 10 9 km. 





A x 


A 2 


^3 


A 4 




A 6 


Ft 


F 2 


T (days) 


17.3 


23 


31 


41 


43.5 


57.5 


27 


30 




0.099 


0.132 


0.177 


0.235 


0.249 


0.329 


0.155 


0.172 


V-MP 


0.13 


0.23 


0.41 


0.72 


0.81 


1.41 


0.31 


0.38 



5.6 Mass moment method 

A more systematic,way to extract information on the masses is to use the method 
of mass moments l3. The mass moments < fi m > are defined as 

< lT >= J dfi e n (fx) ^-fi m . (30) 

< \i m > is related to < r" >= ^ events r n , with t = {vh/te)T, as constructed 
from the observations and which can also be computed as follows 

< r n >= JdN ev e„(/!) T n = Vu TH 7 (m) < M m > , (31) 

with m=(n+ l)/2. For targets in the LMC 7(771) = T(2 - m)H(m) and 

V = 2iV*U s F) rE vh — 2.4 x 10 3 pc 3 N * Us , (32) 

10° star — years 
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(•OO / \ 1—71 

T(2 - m) = J i-^-j f(v T )dv T , (33) 

H{m) = [ (x(l - x)) m H(x)dx . (34) 
Jo 



The efficiency e n (/f) is determined as follows (seeEl) 



( ) = [dKM e(T) r" 



where dN* v (fl) is defined as dN ev in eg. (J10|) with the MACHO mass distribution 
concentrated at a fixed mass /2: dno/dp = uq S(fi — t(T) is the experimental 

detection efficiency. For a more detailed discussion on the efficiency see ref. E 3 ! 

A mass moment < fi m > is thus related to < r" > as given from the measured 
values of T in a microlensing experiment by 



< r" > 
VuTH^ym) 

The mean local density of MACHOs (number per cubic parsec) is < /i° >. The 
average local mass density in MACHOs is < /i 1 > solar masses per cubic parsec. In 
the following we consider only 6 (see Table 3) out of the 8 events observed by the 
MACHO group, in fact the two events we neglect are a binary lensing event and an 
event which is rated as marginal. The mean mass, which we get from the six events 



detected by the MACHO team, is 

0.27 M S . (37) 

(To obtain this result we used the values of r as reported in Table 3. whereas 
r(l)H(l) = 0.0362 and T(2)H(0) = 0.280 as plotted in Fig. 6 of ref. El). When 
taking for the duration T the values corrected for "blending", we get as average 
mass 0.34 M Q . If we include also the two EROS events we get a value of 0.26 M 
for the mean mass (without taking into account blending effects). The resulting 
mass depends on the parameters used to describe the standard halo model. In 
order to check this dependence we varied the parameters within their allowed range 
and found that the average mass changes at most by ± 30%, which shows that the 
result is rather robust. Although the value for the average mass we find with the 
mass moment method is marginally consistent with the result of the MACHO team, 
it definitely favours a lower average MACHO mass. 

One can also consider other models with more general luminous and dark mat 
distributions, e.g. ones with a flattened halo or with anisotropy in velocity space 
in which case the resulting value for the average mass would decrease significantly. 

Another important quantity to be determined is the fraction / of the local dark 
mass density (the latter one given by po) detected in the form of MACHOs, which 
is given by / = M Q /p Q ~ 126 pc 3 < M^-i>. Using the values given by the MACHO 
collaboration for their two years dataE3 (in particular uth = 0.661 corresponding 
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to A > 1.75 and an effective exposure N+tobs of ~ 5 x 10 6 star-years for the observed 
range of the event duration T between ~ 20 - 50 days) we find / ~ 0.54, which 
compares quite well with the corresponding value (/ ~ 0.45 based on the six events 
we consider) calculated by the MACHO group in a different way. The value for / 
is obtained again by assuming a standard spherical halo model. 

Similarly, one can also get information from the events detected so far towards 
the galactic bulge. The mean MACHO mass, which one gets when considering the 
first eleven events detected by OGLE in the galactic bulge, is ~ 0.29M Q Ea. From 
the 40 events discovered during the first year of operation by the MACHO teamE3 
(we considered only the events used by the MACHO team to infer the optical depth 
without the double lens event) we get an average value of 0.16M©. Both values 
are obtained under the assumption that the lenses are located in the disk. For a 
more detailed analysis see ref E3. The lower value inferred from the MACHO data 
is due to the fact that the efficiency for the short duration events (~ some days) 
is substantially higher for the MACHO experiment than for the OGLE one. These 
values for the average mass suggest that the lens are faint stars. 

Once several moments < fi m > are known one can get information on the mass 
distribution dno/dfi. Since at present only few events towards the LMC are at 
disposal the different moments (especially the higher ones) can be determined only 
approximately. Nevertheless, the results obtained so far are already of interest and 
it is clear that in a few years, due also to the new experiments under way (such as 
EROS II, OGLE II and MOA in addition to MACHO), it will be possible to draw 
more firm conclusions. 



6 Formation of dark clusters 

A major problem concerns the formation of MACHOs, as well as the nature of the 
remaining amount of dark matter in the galactic halo. We feel it hard to conceive a 
formation mechanism which transforms with 100% efficiency hydrogen and helium 
gas into MACHOs. Therefore, we expect that also cold clouds (mainly 0LHaI.sb.n1.1M 
be present in the galactic halo. Recently, we have proposed a scenario OHEjIZJO 
in which dark clusters of MACHOs and cold molecular coulds naturally form in the 
halo at galactocentric distances larger than 10-20 kpc, with the relative abundance 
possihliz-jdepending on the distance. Similar scenario have also been considered in 
refs. E30. 

The evolution of the primordial proto globular cluster clouds (which make up 
the proto-galaxy) is expected to be very different in the inner and outer parts of 
the Galaxy, depending on the decreasing ultraviolet flux (UV) from the centre as 
the galactocentric distance R increases. In fact, in the outer halo no substantial 
H.2 depletion should take place, owing to the distance suppression of the UV flux. 
Therefore, the clouds cool and fragment - the process stops when the fragment 
mass becomes ~ 10~ 2 — 10 _1 A/©. In this way dark clusters should form, which 
contain brown dwarfs and also cold i?2 self-gravitating cloud, along with some 
residual diffuse gas (the amount of diffuse gas inside a dark cluster has to be low, 
for otherwise it would have been observed in the radio band). r~ in 

We have also considered several observational tests for our modelEHZi In par- 
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ticular, a signature for the presence of molecular clouds in the galactic halo should 
be a 7-ray flux produced in the scattering of high-energy cosmic-ray protons on H2 ■ 
As a matter of fact, an essential information is the knowledge of the cosmic ray 
flux in the halo. Unfortunately this quantity is unknown and the only available 
information comes from theoretical considerations. Nevertheless, we can make an 
estimate of the expected 7-ray flux and the best chance to detect it is provided 
by observations at high galactic latitude. Accordingly we find a 7-ray flux (for 
E 1 > 100 MeV) $ 7 (90°) ~ / (0.4 - 1.8) x 10~ 5 photons cm" 2 s" 1 sr" 1 (/ stands 
for the fraction of halo dark matter in the form of gas), if the cosmic rays are con- 
fined in the galactic halo, otherwise, if they are confined in the local galaxy group 
Ea $ 7 (90°) ~ / (0.6 - 3) x 10~ 7 photons cm" 2 s" 1 sr _1 . These values should be 
compared with the measured flux by the SAS-II satellite for the diffuse background 
of (0.7 — 2.3) x 10~ 5 photons cm~ 2 s _1 sr^ 1 or the corresponding flux found by 
EGRET of ~ 1.1 x 10~ 5 photons cm~ 2 s _1 sr _1 . Thus, there is at present no con- 
tradiction with observations. Furthermore, an improvement of sensitivity for the 
next generation of 7-ray detectors will allow to clarify the origin of this flux or yield 
more stringent limits on /. 



7 Conclusions 

The mistery of the dark matter is still unsolved, however, thanks to the ongoing 
microlensing and pixcl-lcnsing experiments there is hope that progress on its na- 
ture in the galactic halo can be achieved within the next few years. Substantial 
progress will also be done in the study of the structure of our Galaxy and this es- 
pecially once data from the observations towards the spiral arms will be available. 
Microlensing is also very promising for the discovery of planets. Although being 
a rather young observational technique microlensing has already allowed to make 
substantial progress and the prospects for further contribution to solve important 
astrophysical problems look very bright. 

It has also to be mentioned that it is well plausible that only a fraction of the 
halo dark matter is in form of MACHOs, either brown dwarfs or white dwarfs, in 
which case there is the problem of explaining the nature of the remaining dark 
matter and the formation of the MACHOs. Before invoking the need for new 
particles as galactic dark matter candidates for the remaining fraction, one should 
seriously consider the possibility that it is in the form of cold molecular clouds. A 
scenario this, for which several observational tests have been proposed, thanks to 
which it should be feasible in the near future to either detect or to put stringent 
limits on these clouds. 
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